Casimir forces in multilayer magnetodielectrics with both gain and loss 
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A path-integral approach to the quantization of the electromagnetic field in a linearly amplifying 
magnetodielectric medium is presented. Two continua of inverted harmonic oscillators are used to 
describe the polarizability and magnetizability of the amplifying medium. The causal susceptibilities 
of the amplifying medium, with negative imaginary parts in finite frequency intervals, are identified 
and their relation to microscopic coupling functions are determined. By carefully relating the two- 
point functions of the field theory to the optical Green functions, we calculate the Casimir energy and 
Casimir forces for a multilayer magnetodielectric medium with both gain and loss. We point out the 
essential differences with a purely passive layered medium. For a single layer, we find different bounds 
on the Casimir force for fully amplifying and for lossy media. The force is attractive in both cases, 
also if the medium exhibits negative refraction. From our Lagrangian we also derive by canonical 
quantization the postulates of the phenomenological theory of amplifying magnetodielectrics. 

PACS numbers: 12.20.Ds, 03.70.+k, 42.50.Nn 



I. INTRODUCTION 

The Casimir force is a pure quantum effect that 
can be considered as the macroscopic manifestation of 
the vacuum fluctuations of the electromagnetic fields in 
the presence of boundaries. Originally derived in 1948 
for two ideally conducting or reflecting plates in vac- 
uum |lj , the Casimir force per area Fq j A was found to be 
—hcK^ I (240o?^), an attractive force with characteristic in- 
verse fourth power dependence on the plate separation d. 
A theory for Casimir forces between parallel dielectrics 
was developed by Lifshitz et al. 0, further refined by 
Schwinger et al. [3| , and since then extended to arbitrary 
multilayer dielectrics j^-Q and other geometries. 

In recent years, the Casimir force has become techno- 
logically relevant, with the development of micro- and 
nano-electromechanical systems with small components 
at close proximity of each other. On the one hand, the 
Casimir force can be a major cause of stiction {i.e. mi- 
croscopic components sticking together), friction, or ad- 
hesion, and thus forms a possible obstacle for the oper- 
ation of nanostructured devices. On the other hand, a 
novel class of microelectromechanical systems (MEMS) 
only works because of the Casimir force [1, [13] ■ 

It is therefore both fundamentally interesting and tech- 
nologically relevant to what extent the Casimir force can 
be controlled by changing the electromagnetic environ- 
ment. For most geometries, the Casimir force between 
two media separated by vacuum is an attractive force, 
with a magnitude that becomes appreciable in the submi- 
cron range and rapidly increases in the nanometer range. 
However, already Lifschitz predicted that the Casimir 
force for parallel dielecric layers can be attractive or re- 
pulsive, depending on the relative values of the dielectric 
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constants of the successive layers 0]- The first experi- 
mental confirmation of his prediction came only recently: 
Munday et al. found Casimir forces with different signs 
for suitably chosen interacting materials immersed in a 
fiuid [Tl| . with a measured repulsive interaction which 
was weaker than the attractive counterpart. Further 
measurements of Casimir forces are reported in p^ - [20l |. 

If Casimir forces could be made repulsive, then 
this could eliminate the unwanted phenomenon of stic- 
tion ^l|, enable quantum levitation of objects in a fiuid, 
and lead to new classes of switchable nanoscale devices 
with ultra- low static friction (22l - [25l | . So what options ex- 
ist to make the Casimir force repulsive, besides Lifschitz's 
suggestion? In the first place, metamaterials have been 
proposed to this end. However, loss in metallic substruc- 
tures may intuitively turn a repulsive force into an attrac- 
tive one, which may explain why repulsive Casimir forces 
with metamaterials have not yet been reported [1^ [13] . 
Recently, Zhao et al. showed theoretically that a repul- 
sive Casimir force could be realized with metamaterials 
with strong chirality [2^. Another mechanism to ob- 
tain repulsive Casimir forces is Boyer's Casimir repul- 
sion based on an asymmetric three-layer setup of a non- 
magnetic medium on the one end and a purely rnagnetic 
medium on the other, separated by vacuum [2^, [3^. It 
relies on the nontrivial possibility of developing new ar- 
tificial negative- index metamaterials [27]. 

A third option to obtain repulsive Casimir forces, and 
the one considered in this paper, is the use of media with 
optical gain [Sllll^- Indeed, the main aim of the present 
article is to determine Casimir energy and forces in am- 
plifying magnetodielectrics. With amplifying medium, 
we mean a medium for which the imaginary part of the 
electric or the magnetic susceptibility becomes negative 
for one or more frequency intervals (Im[e(aj)] = ei(w) < 
or /ii(w) < 0, in contrast to lossy systems for which both 
ei[(jj) and /ii(w) are always positive. We allow for gain 
not only in the electric but also in the magnetic response. 
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thereby treating electric and magnetic fields on an equal 
footing in our theory. 

Since it is a quantum mechanical effect, the calcula- 
tion of the Casimir force for media with gain requires 
a consistent procedure for quantization of the electro- 
magnetic field in presence of amplifying medium. This 
can be done with the concept of inverted quantum har- 
monic oscillators that was introduced by Glauber [33| . 
Generally, a rigorous quantization procedure would re- 
quire a Lagrangian and Hamiltonian formulation of the 
theory, followed by the standard canonical quantization 
rules. Currently, a consistent phenomenological approach 
exists to macroscopic quantum electrodynamics in pres- 
ence of amplifying media, but no canonical formulation 
is attempted |34| - l38j . Below we derive the postulates of 
this phenomenological theory from a canonical quantum 
theory, where continua of inverted oscillators are used to 
describe linearly amplifying media. 

The Casimir force in the presence of amplifying materi- 
als, with or without negative index, has only recently be- 
gun to be explored. Leonhardt and Philbin calculated the 
effects of an amplifying dielectric on the Casimir force, 
based on the assumption that the well-known Lifshitz for- 
mula for the Casimir force is applicable without change 
to amplifying media [32i] . and found that the Casimir 
force in the presence of this medium is repulsive. We will 
address this issue here as well. 

Sambale et al. [H, H^] apply the phenomenological 
quantization of the electromagnetic field for amplifying 
magnetodielectric media to calculate the Casimir and 
Casimir-Polder forces. They find that the Casimir- 
Polder force on a weakly polarisable plate of excited gas 
atoms is attractive at short distances from a mirror, and 
oscillating behavior between attraction and repulsion for 
larger plate-mirror separations. However, at a more tech- 
nical level there is some controversy about the applica- 
bility of the Minkowski stress tensor and an alternative, 
Lorentz force-based tensor was proposed [4]| . This ap- 
proach was adopted in Refs. 39, 1^ to calculate the 
Casimir(-Polder) forces, but disputed by Pitaevskii [l^l 
and Brevik 

Our approach to calculate Casimir forces is different. 
We first develop a path-integral method for the quanti- 
zation of the electromagnetic field in linearly amplifying 
magnetodielectrics. We benefit from and generalize re- 
cent results obtained with path integrals. In Ref. [3], 
Li and Kardar developed a path-integral approach for 
computing fluctuation-induced forces between manifolds 
immersed in a correlated fluid. Golestanian and Kardar 
extended this formalism to arbitrary but small deforma- 
tions of the boundaries and focused on the mechanical 
response of the vacuum [4^. Emig and colleagues also 
used the path-integral formalism to obtain the normal 
and lateral Casimir forces between two sinusoidally cor- 
rugated perfectly conducting surfaces [i^ . Recently, one 
of the authors extended this formalism to calculate the 
Casimir force between two perfectly conducting plates 
immersed in a magnetodielectric medium ,47j . Here we 



extend this quantization scheme to arbitrary multilayer 
amplifying media. 

The structure of the paper is as follows: In Sec. [Hi 
we propose a Lagrangian for the electromagnetic field 
in an amplifying magnetodielectric medium and derive 
the generating function, which is used in Sec. Illll for the 
path integral quantization. Causal electric and magnetic 
susceptibilities of the amplifying medium are obtained, 
both in the frequency interval(s) with gain and in the 
remaining lossy regions. We calculate the Green tensor 
and Casimir forces for a multilayer amplifying magne- 
todielectric medium in Sec. IIVI and corresponding nu- 
merical results are presented in Sec. [V] In Sec. I VII we 
derive from our Lagrangian a canonical theory that sup- 
ports the known phenomenological approach to quanti- 
zation in amplifying magnetodielectrics. We conclude in 
Sec. I VIII Further details of our calculations are given in 
two appendices. 



II. FIELD QUANTIZATION 

The quantum electrodynamics of a linearly damped 
magnetodielectric medium can be described by modeling 
the medium as two independent reservoirs that interact 
with the electromagnetic field. Each reservoir contains 
a continuum of three-dimensional harmonic oscillators 
that describe the polarizability and magnetizability of the 
medium [48l - l53 |. We assume the medium is linearly am- 
plifying, at least in one or more finite frequency windows 
where ei{uj) < and/or /ii(aj) < 0. Despite the striking 
practical differences between amplifying and lossy media, 
their theoretical descriptions turn out to be quite simi- 
lar. We adopt Glauber's inverted oscillator to model the 
quantum amplifier [33,] and use continua of inverted oscil- 
lators to describe gain instead of loss in the polarizability 
and magnetizability of the medium. 

We introduce our model for optical media with both 
gain and loss by first specifying its Lagrangian density in 
real space, 

^ — .CeM 4- £e + + ^int, (1) 

where the electromagnetic part £em has the standard 
form £em = 5eoE^(x, i) — 2^B^(x, t). There is gauge 
freedom to write the electric field E = —dA/dt — Vcj) and 
the magnetic field B = V x A in terms of the scalar and 
vector potentials (j) and A. For convenience we choose the 
Weyl gauge in which the scalar potential vanishes, which 
allows us to write E and B in terms of only the vector 
potential. The amplifying magnetodielectric medium is 
modeled with frequency continua of independent vector 
fields X^^ (x, t) and Y^^ (x, , in terms of which we will de- 
scribe the linear electric and magnetic polarization of the 
medium. Therefore the material part of the Lagrangian 
density describing the amplifying medium can now be 
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written as 



r — 



-X5(x,i)--a.^X^(x,t) 



sgn[ei(w)], 
(2a) 

sgn[/ii(w)]. 

(2b) 



In these Lagrangian densities, for frequencies with pos- 
itive signs the medium is fossy, and otherwise it is am- 
pUfying as modeled with inverted osciUators due to the 
minus sign. As we wiU see in Sec. IVI[ this modifica- 
tion leads to field operators that satisfy Maxwell's equa- 
tions and whose positive-frequency components are as- 
sociated with both annihilation and creation operators 
in the case of amplifying media, in accordance with the 
previous works (26-38] . 

We define the polarization and magnetization fields of 
the medium as 



P(x,t) = 
M(x,t) 



dw/(x,w)X^(x,t), (3a) 

O 

dwg(x, w)Y„(x,t), (3b) 



and assume a linear coupling of the electromagnetic field 
with these fields, 

£i„t(A, P, M) = A(x, t) ■ P(x, t) + V X A(x, t) ■ M(x, t). 

(4) 

The /(x, w) and g{x,uj) in Eq. ^ are the real-valued 
scalar coupling functions of the inhomogeneous medium 
and the electromagnetic field. We have implicity assumed 
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that the medium is isotropic by taking scalar coupling 
functions. Anisotropy could be included by making them 
tensors, but this is not pursued here. 

With the Lagrangian in a suitable form, we can now de- 
fine a generating function for our path-integral quantiza- 
tion. For a field theory with only a single scalar canonical 
field if, the generating functional (or partition function) 
has the form [551 



Z[J]= / V[ifi]exp 



d'^x[C{ip{x)) + J{x)ip{x)] 



(5) 

where a; G is a spacetime coordinate, and J is the 
auxiliary source field associated with the scalar field ip. 
In our case we have several interacting canonical fields in 
our total Lagrangian (jlj , so we need to generalize Eq. ([5]) . 
We first calculate the partition function Zq for the free 
fields, i.e. neglecting their interactions: 



Zo[J 



EM 



X exp 



1^ J ^'^^ i^EM 



V[A]V[X^]V[Y^] 

£c + An + Jo - A 



did Jp 



(6) 



in terms of the auxiliary source vector field Jem foi' 
the electromagnetic field, and the frequency continua 
of source fields 3c,uj and Jni,tj associated with the elec- 
tric and magnetic polarization fields, respectively. We 
can write Zq in a more convenient form by employing 
the four-dimensional version of Gauss's theorem for the 
vector potential, thereby replacing £em in Eq. © by 
[— ^ A • (V X V X A) — EqA ■ dfA] , and by using inte- 
gration by parts for the polarization fields X^^ and Y^^. 
This gives 



■^o[Jem, Jc 



2?[A]2?[X„p[Y„]exp^-- 



d^^ 



2 A • (a^o" 



X Vx 



1 



duj ( Jc,tj(x) • X^(x) 4 

uj^)Y^{x)ag-a[fii{uj)] 



^Ux)-Y^{x) + -XUx)-{d, 



- eodf)A- 3em{x) ■ A{x) 



a;2)X„(a;)sgn[ei(cj)] 



The above partition function is Gaussian since the inte- 
grand is quadratic in terms of the fields. Therefore the 



J 



functional integration can be performed exactly and the 
result is .^o[Jem, Jc,tj, Jin,tj] = exp{i£o/2fi}, in terms of 
the Lagrangian for non-interacting fields 



£0 d*x J d^x'JEuix) ■ G'^l{x~^x')-3EMix') + Jd^xJdtJ dt' J dw 



{Jc,c^(x, t) ■ Gc,c^(t - t') ■ Jo,c^(x, t')sgn[ei{uj)] + Jm.Lui^, t) ■ G,n,c^(i - t') ■ Jm,i^(x, i')sgn[/ii (w)] } . (7) 



Here, the space component of the space-time x G is indicated in bold by x G M'^ and the time component 
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by f G K. The Green tensor 
electromagnetic field satisfies 



,(0) 

'em 



{x — x') for the free 



V X V X 



1^ 

C2 9t2 



G^°' (x^ 



■x')^^ioS\x-x'), (8) 



describing the propagation of light in time and in free 
space, whereas the Green tensors Ge,u{t — t') = G^^^{t — 
t') for the non-interacting electric and magnetic polariza- 
tion fields describe propagation only in time, 



q2 

— + I G,/^ - t') = 5{t - t')U 



(9) 



where I3 denotes the spatial miit tensor. This indi- 
cates that the only way to transport energy in the in- 
teracting system is via the electromagnetic field. For the 
same reason, the polarization and magnetization fields 
in the absence of the electromagnetic field do not lead 
to a Casimir force [13] . The retarded solution of Eq. © 
in Fourier space is Ge/i„,tj(w') = l3/[a;^ — [to' + iO+)^]. 
It is convenient to define the source fields Jp^M for the 
electric and magnetic polarization fields as linear com- 
binations of the corresponding frequency continua Je,m, 
namely Jp(x) = dw /(x, aj)Jo,(.j(x, cj) and Jm(x) = 
dti;5(x, w)Jni^ij(x, w). With these, the generating 
functional for the interacting fields can be written in 
terms of the free generating functional as [ssj 

Z[Jem,Jp,Jm] =^(r'[0,0,0] 

/- [ A'^xC- ^ ft 5 h 6 

\h J \ i dJEuix) ' i d3p{x) ' i dJuix) 



xZo[Jem, Jp, J 



MJ 



(10) 



where Zo[Jem, Jp, Jm] is the fi:ee-space partition func- 
tion and £int is given in Eq. ^ and Z^j"^ [0,0,0] is nor- 
malization factor. The exponential in this functional is to 



be understood as a power series in the coupling functions, 
that is by perturbation theory. By using the specific form 
Eq. (|4]) for the interaction, we obtain 



Z[Jem,Jp,Jm] = Ztr'[0,0,0]^- 



n=0 



Ah 



d X 



d 



(5 Jem (a;) dtS3p{x) 
xZoJJem, Jp, Jm]- 



V x 



(11) 



Hereby we determined as one of our main results the par- 
tition function Z for the interacting fields that describe 
a magnetodielectric medium with both linear gain and 
loss. 



III. GREEN TENSORS AND 
SUSCEPTIBILITIES OF THE AMPLIFYING 
MEDIUM 

As we will see below, the Casimir force can be com- 
puted in terms of the electromagnetic Green tensor of the 
medium. In our zero-temperature field theory, the Green 
tensors (or propagators) are vacuum expectation values 
of time-ordered products of field operators, which can be 
computed as functional derivatives of the partition func- 
tion (see also [U 



>em 



(x, x') = —ih 



(^^ZJJem, Jp, Jm] I 
(5JEM(a;)^JEM(a;')'''™= 



.o'(12) 



in terms of spacetime coordinates x,x' . The medium de- 
scribed by our Lagrangian ([l} in general is not transla- 
tionally invariant, but it is stationary, and consequently 
GEM(a;,a;') = Gem(x, x',i — t'). After evaluating the 
functional derivatives of Eq. ([T^ we obtain a Dyson equa- 
tion for the Green tensor that after time-Fourier trans- 
formation becomes 



GEM(x,x',a;) = g'^1^{x - x.' , uj) 



a;^y"dxi[G^°]^(x-xi,a;) • y dtj'{sgn(ei(a;'))/^(xi, iu')Go,c^'(a;)} • Gem(xi, x', w)] 

y dxi[G|5%(x- xi,w) X Vi • y di:j'{sgn(^i(w'))5^(xi,iLj')Gni,a;'('^)} • Vi x Gem(xi,x',w)]. (13) 



This long equation for the Green tensor can be brought 
into a more familiar form by applying the differential 
operator {fiQ^V x V x — w^eols) to both sides, giving 

V X [/Lt"^(x,a;)V x Gem(x, x', w)] 

-^^^^^^GEM(x,x',c.)=/io<53(x-x')l3. (14) 

Here we defined the electric permittivity e(x, w) — 
1 + Xc(x, w) and the inverse magnetic permeability 



/i ^(x, w) = 1 — Xm(x, w) of the amplifying magnetodi- 
electric via 

Xo(x,a;) = — / duj — — . , (15a) 



Xm(x, w) = flQ 



duj' 
duj' 



So Jo (w-|-i0+)2 ' 

,.g2(x,w')sgn[Aii(a;')] 



(15b) 



/o a;'2 - (cj + i0+)2 

With these definitions and the requirement on the 
coupling functions that /^(x, — oj*) = /^(x, a;) and 
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(/^(x, —a;*) = (/^(x, uj), the e and fi are complex functions 
of frequency which satisfy Kramers-Kronig relations |66j] 
and have the properties of the response functions i.e, 
e(x, —u!*) = e*(x, Lu), and analogously for /i. Purely lossy 
media would have the further properties £i(x, w) > 
and /ii(x, oj) > 0, but here we have a model that can 
describe amplification in some frequency interval(s) as 
well, for which ei(x, cj) < and/or /xi(x, w) < 0. The 
functions Xc,m(x, w) have no poles in the upper-half fre- 
quency plane and tend to zero as — > oo, so that in the 
time domain, the electric and magnetic susceptibilities 
Xe,m(x, t) corresponding to Eq. (|15p become proportional 
to the step function <d{t). This is as it should be, since 
either with gain or loss, the response should be causal. 

It is important to stress the generality of our model: 
if we are given definite functions for the electric per- 
mittivity £(x, w) and magnetic permeability /i(x, w) of 
the gain medium, then we can invert the relations (jlSp 
to find the corresponding coupling functions /(x, w) = 

^2weo|ei(x,a;)|/7r and g{x,uj) = ^2a;|/ij"^(x,a;)|/7r^o, 
where the modulus signs ensure that the coupling func- 
tions are real-valued both for lossy and for amplifying 
media. A similar general theory, albeit for purely lossy 
media, can be found in Refs. HI, Specific choices 
for the optical functions e and fj, will be made for our 
numerical investigations in Sec. IVl below. 

We will make some further consistency checks on 
the path-integral quantization for amplifying dielectrics. 
Similar checks for lossy media were performed in 
Ref. pf?! . Recall that the defining equation (fH|) for elec- 
tromagnetic Green tensor is found by functional differ- 
entiation of the partition function. This equation en- 
ables the identification of the dielectric functions e{u!) 
and ^(w) for the amplifying medium. Analogously, we 
can find Green tensors for the material fields in our the- 
ory, as well as correlations functions of mixed type. An 
example of the latter type is 



Gem,p(x,x', w) 



-ih 



(5^Z[Jem, Jp, Jm] 



6JEMix)S3p{x') IJem,Jp,Jm-0 

icjEo [e(x,t^) - 1] GEM(x,x',a;). (16) 



Analogously we find Gem,m(x, x', cj) — 
^0 ^ [l - ^~^(x, w)] V X Gem(x,x',iu), with £{uj) 
and fi{i-u) as previously defined in Eqs. (jlSp . This 
shows that the definition of these response functions for 
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amplifying media can be made uniquely and consistently 
in the path- integral quantization method. In Sec. IVII we 
will derive the equivalent canonical quantization theory 
for amplifying dielectrics, with the same Lagrangian ([T]) 
as a starting point. 

IV. CASIMIR FORCE FOR AMPLIFYING 
MULTILAYER MEDIA 

A. Derivation of Fc in path-integral formalism 

Here we calculate the Casimir force for two parallel 
perfectly conducting plates that are separated by a mul- 
tilayer linearly amplifying medium of total width d. Of 
course perfect conductors do not exist, and the assump- 
tion of linear amplification in reality will break down in 
an amplifying medium without round-trip losses. Still, 
an important advantage of our model is that we consider 
causal optical response functions e{uj) and satisfy- 
ing the Kramers-Kronig relations, as it should for any 
medium that respects causality, amplifying or not. 

For the single homogeneous amplifying medium, the 
Casimir force between two plates can be computed as 
the spatial derivative of the effective action 



d 



(17) 



where the effective action is proportional to the logarithm 
of the partition function. 



S'cff(d) = h\nZ[d]. 



(18) 



Now since for planar structures there are independent 
TE and TM solutions of Maxwell's equations, the total 
partition function is the product of ZxEic^] and .ZtmI^] 
partition functions, so that the effective action intuitively 
becomes the sum of TE and TM contributions 



Scs[d] = ^(IuZtmM] + IuZteM]), 



(19) 



and likewise for the Casimir force. The details of the cal- 
culation of the partition functions are left to AopendixEl 
and yield 



1 



^TE,TM 



det Fte.tm {x,y,zi,Z2)' 



(20) 



where 



rTM(a;, y, ^1,^2) 
TTE{x,y,zi,z2) 



Gtm{x - y,zi,zi) Gtm{x - y,Z2,zi) 
Gtm{x - y,zi,Z2) GTuix - y,Z2,Z2) 

-d'^GTEix - y, zi,zi) -d'^GTEix - y, Z2, zi) 
-d'^GTEix - y, zi,z2) -d'^GTEix - y, Z2, Z2) 



(21a) 
(21b) 



r 



We see from Eqs. (|17ll2ip how the Casimir force is ex- pressed in terms of Green tensors ^te,tm, which can be 
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obtained from the Green tensor Eq. (fH]) by applying a 
Wick rotation. The explicit form of the Green tensor 
for planar multilayer dielectric structures were obtained 
in Ref. [56| . The details of the calculation of the Green 
tensors for the more general situation of amplifying mag- 
nctodielectric multilayer media are summarized in Ap- 
pendix |B] 

Before considering multilayer media in more detail be- 
low, we focus on Casimir force in the presence of a sin- 
gle homogeneous amplifying layer. This is the geome- 
try originally studied by Casimir, but now with the vac- 
uum between the conductors replaced by the amplify- 
ing medium. For this simple geometry, the Dirichlet and 
Neumann boundary conditions are formally the same and 
lead to the same result, so that TE and TM waves each 
account for half of the Casimir force. 

For this case, the Green function ([Ti)) in 2D 
Fourier space and after Wick rotation can be writ- 
ten as gEM( q, i^, z') = fiofi{iuj) e '/(2Q), where 
Q(q, iw) = y/q'^ + uj'^e{iuj)^{iuj)/c'^ (see Appendix [B]) . 

There is no ambiguity how this square root is to be 
taken, since e{iuj) and fi{iuj) are both real-valued func- 
tions of (real) w. Since both e(w) and //(w) have no 
odd-order zeroes in the upper-half frequency plane and 
tend to unity in the limit of going to infinity, it fol- 
lows that both e and fi assume positive real values on 
the positive imaginary frequency axis [57i] . Since conse- 
quently e(a;)/i(cj) does not have any poles or odd-order 
zeros in the upper half-plane of frequency, then n{uj) for 
Im(cLi) > is defined as the analytic branch of e{u!)fi{uj) 
that tends to +1 as |a;| goes to infinity. (Otherwise, 
■y/ e{uj)fi{uj) would not be an analytic function there and 
corresponds to materials with so-called absolute instabil- 
ities [H, [5^.) So we find that both n{uj) and the Green 
tensor 5EM(q, w, z, z') are analytic in the whole upper 
complex- frequency plane, as it should |6Qi] . By substitut- 
ing the expression for the Green function Gem{ci, iw, z, z') 
into Eq. (IT51) . the Casimir force per unit area for a ho- 
mogeneous amplifying medium becomes 



Fc=-h 



duj(fq 2Q(q,iw) 
(27r)3 e2S(q,ic>;)<i _ 1 



n{iuj)uj ] 



do; 



(22) 



The last identity in (1^^ follows from a partial integra- 
tion over uj an d up on calculating the w-derivative of the 
integral over q [6l| . For the empty cavity, the integrals in 
Eq. ((22|) can be evaluated exactly, giving the well-known 
resuh Fc = -hcn^ /{2i0d'^). 

As one of our main results, we find that the Casimir 
force in the presence of an amplifying medium (j22p has 
the same form as for a purely attenuating medium (47| . 
provided that e{uj) 11(00) for amplifying medium does not 
have any poles or odd-order zeros in the upper half plane. 
Thus the actual value of the force for lossy and amplifying 
media can only follow from the different forms of s{uj) 
and /i(a;) in both cases. The key difference was already 



stated in the Introduction, namely that amplifying media 
have one or more frequency intervals with ei{uj) < or 
/ii(a;) < or both, for real frequencies a; > 0, whereas 
lossy media always have £i(w), fJ-i{io) > for positive real 
frequencies. 



B. Analytical results: bounds on the Casimir force 

Realistic amplifying media are amplifying in one or 
more frequency intervals and lossy elsewhere. Since the 
Casimir force is obtained as an integral over all fre- 
quencies, it may well be that the lossy part dominates 
the total Casimir force. We will study these issues nu- 
merically in Sec. IV Bl To understand the effect of am- 
plification on the Casimir force, we will first make the 
further assumption that the medium is fully amplifying, 
by which we mean that ei(w), /xi(a;) < not only for 
some but for all positive frequencies. We make this ad- 
mittedly unrealistic assumption to single out the effect of 
linear amplification on the Casimir force. This will give 
us some insight, and after that in Sec. |V] we will relax 
the assumption of full amplification. 

Our aim here is to give a bound for the Casimir force in 
the presence of a fully amplifying medium, similar to the 
bounds obtained in Ref. 62| for passive dielectric (i.e. 
nonmagnetic) mirrors, where from causality considera- 
tions it follows that the Casimir force on dielectric slabs 
is always attractive, but less so than between two ideal 
mirrors separated by vacuum. 

The force ([^^ depends on the refractive index n{iuj) — 



^ e(iw)/i(iaj), in terms of the susceptibilities e and /i that 
both tend to unity for high frequencies. Causality implies 
the identity [57] 



e[ibj) 



(23) 



and an analogous identity holds for /i(w). It follows and is 
known that for passive systems, with £1(0;) and /ii(w) al- 
ways positive on the positive real frequency axis, e(iw) 
and /^(iw) decrease monotonically from a finite value 
^static, ^static > 1 (or +00 for the electric response of 
metals) at w = down to unity for a; — >■ oo [stI ]. 

We will instead apply the identity (^5)) to amplifying 
media, for which it also holds as long as they are de- 
scribed by causal response functions. It follows that for 
fully amplifying systems, with £i(w) < and /ii(w) < 
on the whole positive real frequency axis, siiuj) and /^(iw) 
from a finite value < 1 at a; = increase monotonically 
towards unity for w — > oo. We can say even more, us- 
ing that causal response functions e(w) and /i(w) have 
no zeroes in the upper-half frequency plane [57]: for fully 
amplifying media we find that £(iuj) and /i(icj) increase 
monotonically from finite values < Estatic, Mstatic ^ 1 
towards unity for cj — )■ oo. 

How are these results related to the Casimir force? 
We have just found that for fully amplifying media. 
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n{iu!) = e{iuj)^{iuj) increases monotonically, assuming 
values between ngtatic ^ and 1. Therefore, there is a 
1-to-l mapping 

a; -!-> s = ujn{iLj), w, s € [0, oo). (24) 

This allows us to rewrite the Casimir force for fully am- 
plifying media of Eq. (I22p in terms of the new variable s 
as 

Fc = — ^/ ds , ,/ \ . 25 

Notice that for free space we have n(iw) = 1 and 
hence w(s) = s, which immediately gives the well-known 
Casimir force for vacuum, F^'^ — — ftc7r^/(240d'*). But 
also for the general case we can say more: since 
< fistatic < n{iuj) < 1 for fully amplifying systems, we 
can invert the relation (j24p and find s < uj{s) < s/?istatic 
for all s. Combining this with Eq. (j25p immediately gives 
for all separations d the inequalities 



ffl(tE>) 
















4 > 


< > 




< ^ 



z 



FIG. 1. A planar multilayer magnetodielectric medium sand- 
wiched between two perfect mirrors, with parallel planar 
layers labeled hy I = 1, 2, . . . , A. The coordinate system is 
chosen such that the layers are perpendicular to the z-axis. 
Layer j has thickness dj. Each layer is assumed to be ho- 
mogeneous, isotropic and of infinite transverse size such that 
e(x,a;) = ej{ijj) and /i(x,aj) — ^J.j{uo) for x in layer j. The 
conductors have coordinates zi and Z2, and their separation 
d equals dj . 



rpvac 

^ ^tullamp ^ ^vac 

^static 



hcTT 

' 240^4 ■ 



(26) 



In other words, the Casimir force on two ideal conductors 
separated by a fully amplifying medium of width d, is al- 
ways more attractive than if the medium were vacuum. 
But it is not more attractive than by a factor l/ristatic- In 
particular, we find no sign change in the Casimir force (no 
Casimir repulsion) on ideal conductors separated by a ho- 
mogeneous fully amplifying magnetodielectric medium. 
These bounds also hold for fully amplifying media that 
for some frequencies exhibit negative refraction, as nu- 
merical examples in Sec. |V] will illustrate. 

The bound (j26p holds more generally for magnetodi- 
electric media for which niiuj) increases and < n{iuj) < 
1. For example, if /i(a;) is purely lossy and e{uj) describes 
full amplification, then the product of the monotonically 
decreasing fi{iuj) and the monotonically increasing e{iuj) 
may still be a monotonically increasing function between 
and 1. 

Similarly, for passive media the causal response func- 
tions s{uj) and fJ-{uj) have no zeroes in the upper-half 
frequency plane and e(iw) and ^(iw) decrease monoton- 
ically towards unity for cj — ^ cc [s^l- Consequently 
s/ristatic < w(s) < s for all s, and we find the follow- 
ing inequalities for the Casimir force 



(27) 



So the Casimir force on two ideal conductors separated 
by a lossy medium is always attractive, and less attrac- 
tive than in vacuum, but the force is not reduced by a 
factor larger than l/ristatic- These bounds also hold for 
passive media that for some frequencies exhibit negative 
refraction, as numerical examples in Sec. |V] will illustrate. 



C. Casimir forces in amplifying multilayer 
magnetodielectric media 

Here we generalize our previous results for homoge- 
neous media to planar multilayer geometries, with N 
parallel planar layers labeled by ^ = 1, 2, . . . , of thick- 
nesses di, as depicted in Fig. [T] Each layer is assumed 
to be homogeneous, isotropic and of infinite transverse 
size. As is well known for such a planar multilayer geome- 
try, the electromagnetic field can be completely expanded 
into independent transverse-electric (TE) and transverse 
magnetic (TM) fields, that satisfy the same scalar wave 
equation, but differ in their boundary conditions. The 
Green tensor, that via Eq. (PT|) determines the Casimir 
force, can also be separated into TE and TM parts. For 
multilayer dielectric media, the Green tensor was ob- 
tained by Tomas [s^ , essentially using a transfer matrix 
approach, and a generalization to lossy magnetodielec- 
tric can be found in [63l |. For our purposes we need a 
further generalization, namely the Green tensor for am- 
plifying magnetodielectric multilayer media, and in Ap- 
pendix |B] we give a brief derivation and the final result. 
In general, one finds A^^ expressions for the Green tensor 
Q^^''^^{ci,iuj, z, z'), depending on which of the A^ layers 
the two coordinates z and z' are in, but for the Casimir 
force on the two ideal conductors, Eq. (PT|) shows that 
we fortunately only need four of those terms, namely the 
ones for which z and z' both coincide with one of the 
coordinates zi (boundary of layer 1) and Z2 (boundary 
of layer A^) of the ideal conductors. 

Here we will focus on the Casimir force on two ideal 
conductors separated by three slabs of matter with linear 
gain and loss, which are spatially homogeneous in the 
layers 1, 2, and 3 (see Fig. [T]). Following the method 
outlined in Appendix|Bl whether we find the four relevant 
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expressions for the Green functions, both for TE and TM 
polarizations. We give two of them: For z, z' both in layer 
1, we find 



/-"TETM/ • i\ I- \ T-TE,TM 



2Qi 

while for z in layer 1 and z' in layer 3 we obtain 



,(28) 



/-■TETM/ ■ i\ I- \ T-TE,TM 



,TE,TM 



where "ij^ 



1 and 



2Qi 



(29) 



Also, the reflection 



1 , ■ • i-i J TE.TM , ,TE,TM , 

and transmission amplitudes r^_^_ and t^^^ can be 



calculated with the recursive relations (|B6p and (|B7p . 
The other two Green functions can be found analogously, 
for both polarizations. From Eq. (|19p we then flnd the 
effective action for the amplifying three-layer magnetodi- 
electric medium. 



•Soft 



dcjd q 



^[ln{(l + r™e-2S...)(l + ,^^^ 
(TM ^ TE)]. 



(30) 



One can check that in the absence of ampliflcation, our 
exp ression (1301) tends to a known result for lossy media [2i, 
|64| . In the limit where the two perfect conductors are 
brought to infinity (i.e., di^d^ — > oo), Eq. ([5D|) reduces 
to 



'S'cff = ft 



dwd q 



;27r)3 
(TM^ TE)]. 



[ln{(l-r-r-r)e 



,TM TMn -2Q2d2 



} 



(31) 



This is the generalized Lifshitz formula for the Casimir 
energy density for amplifying media, in the specific three- 
layer geometry where two semi-infinite media with per- 
mittivities and permeabilities ei, fj,i and £3, fi^ are sep- 
arated by a medium of permittivity £2 and permeability 

The most realistic special case of Eq. (pij) is the one 
where the two semi-infinite media are lossy rather than 
amplifying, and gain occurs for some frequencies in the 
middle layer. 



V. NUMERICAL RESULTS 
A. Numerical results for fully amplifying media 

In our numerical investigations, we first choose homo- 
geneous single-resonance Lorentz-oscillator models [65] 
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FIG. 2. Frequency dispersion of the real (upper panel) and 
imaginary (lower panel) parts of the electric permittivity and 
the magnetic permeability, with (real) frequencies on a log 
scale. The media are assumed to be single-resonance media 
with resonance frequency ujq and with £{io) = fJ.{ijj) of the 
form described by Eq. (|32p . Parameters are: for the loss 
material (dotted curve): ajp(o,m)/i^o = 0.9, 7{o,m)/i^o = 0.001, 
and for the gain materials the parameters are similar to the 
loss material but with a;p(e,in)/'^o = 0.9 (dotdashed curves), 
and tJp(e,m)/'^o = 0.5 (dashed curves). Note that there are 
frequency regions where refractive index become negative. 



both for the electric permittivity £(w) and for the mag- 
netic permeability n{oj). For simplicity we also assume 
that the electric and magnetic responses are the same. 



e{(jj) = fi{uj) = 1 



9 9 • ' 



(32) 



wher 



is the coupling frequency, ujq the transverse 



resonance frequency, and 7 the amplification parameter. 
The minus sign in front of the second term in Eq. (j32p 
accounts for optical gain that arises from population in- 
version in the medium, and it differs from the usual pos- 
itive sign for passive systems, consisting for example of 
two-level systems in their ground states. The model ([5^ 
is an example of a fully amplifying medium (introduced 
in Sec. lIVBp . since ei{uj),ni{uj) < for all positive fre- 
quencies, as depicted in Fig. [2] 

In Fig. [21 the upper panel depicts the real parts 
£r(iu;),^r(w). The important thing to notice is that 
for frequencies {ujq — ujpY^'^ < u) < ojo both £R(a;) and 
fi-[i{u)) are negative for amplifying media, whereas for pas- 
sive media they are negative in the region < ^ 
(wg -I- Wp)-^/^. So our model with s{uj) = fj,{u}) describes 
negative refraction, and in particular the perfect-lens sit- 
uation n{uj) = — 1 occurs, or at least Re[n(a;)] = — 1. 
The lower panel of Fig. [2] depicts the different signs of 
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FIG. 3. Casimir force Fc per unit area on two perfectly con- 
ducting plates, as a function of their separation, for different 
fully amplifying and lossy dispersive materials between the 
plates. The parameters are identical to those used in Fig. [J] 
The solid curve corresponds to vacuum. The force is given in 
units of hio^lc? and the separation in units of c/a;o. The inset 
shows e = at imaginary frequencies. 
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FIG. 4. The real and imaginary parts of the electric per- 
mittivity and the magnetic permeability in real frequency for 
the mixed-type gain material with electromagnetic parame- 
ters described by Eq. (|33j ; cjq/^^o = 1-5, ijJp/i^o = 0.85, 
uj'Jco^ = 1.8, iV^o = 0.01, = 10. 



£i(aj), for lossy and for fully amplifying media. The 

requirement that e{uj) and have no simple zeroes 

in the upper-half frequency plane does not restrict any 
parameters of Lorentz-oscillator models describing loss, 
but for the fully amplifying model (j32[) we must require 
Wp < Wo [1^. This requirement is usually met, since for 
natural materials the permeability of the medium typi- 
cally equals unity and 70 ^ Wp o ^ woe [ll] . 

In Figure [3] we compare the Casimir force on two per- 
fect planar conductors as described by Eq. (l22t for lossy 
and for fully amplifying homogeneous media, as well as 
for vacuum, all as a function of plate separation d. The 
inset[3]shows the real- valued e{iuj), which is indeed mono- 
tonically decreasing as discussed in Sec. IIVBI here from 
^static = l+Wp/wQ down to unity for lossy media. For am- 
plifying media e(ia;) is indeed monotonically increasing, 
from Estatic = ^ ^ ^p/^o towards unity. For some ampli- 
fying negative-index geo metries, the Casimir force was 
found to be repulsive [S^], but Fig. [3] illustrates our find- 
ing of Sec. lIVBl for a single homogeneous layer: attractive 
Casimir forces for all plate separations for all homoge- 
neous fully amplifying Kramers-Kronig media, including 
the media where for some frequencies there is negative 
refraction. More specifically, it is easily verified that the 
bounds of Eq. (j26|) hold in Fig. [3j the attraction for fully 
amplifying media is always stronger than for plates sep- 
arated by vacuum, but weaker than l/rigtatic times the 
free-space value. Likewise, the analogous bounds on Fc 
in Eq. ((27)) for passive media are also seen to hold in 
Fig. El 



B. Numerical results for homogeneous medium 
with both gain and loss 

The fully amplifying media studied in Sees. IIV Bl and 
IV Al do not occur in nature, but give insight in the ef- 
fect of amplification on Casimir forces. Real amplifying 
media are typically amplifying in a limited frequency in- 
terval, and lossy elsewhere. Therefore we will now study 
the effect of gain in a limited frequency interval on the 
Casimir force, but we still assume = ^J■i^) tor sim- 
plicity. We modify the single-resonance model ([5^ for 
fully amplifying media by adding a loss term. 



(C^f )2 - - ic^7g " (C^l )2 _ _ 

By our choice of parameters, we describe a medium with 
gain around Uq and loss elsewhere, as seen in Fig. 2] 

In Figure [5l we show the corresponding Casimir force 
for this medium, as calculated with Eq. (1^^ . More pre- 
cisely, the figure depicts the difference of the Casimir 
force with respect to the free-space value, and such that 
a positive value corresponds to a more strongly attrac- 
tive Casimir force than for free space. For lossy me- 
dia, we know that this will result in a curve entirely be- 
low the horizontal axis, while for fully amplifying media 
only positive curves would result, as seen in the figure. 
For the medium ([55)1 with both gain and loss and with 
£{uj) = it follows from the inset of Fig.[S]that n{\uj) 

is not monotonically increasing or decreasing, but shows 
more complex behavior: for low frequencies, the curve de- 
creases as for a purely lossy medium. After going through 
a minimum, the curve increases for a while as fully am- 
plifying media would do monotonically, and then finally 
it decreases again. 
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FIG. 5. The difference of the Casimir force with respect to 
the vacuum value, for several dispersive media: lossy, fully 
amplifying (Eq. (|32)l 'l. and mixed- type (Eq. (|33)l '). Scaling of 
the force and of the distance, and parameters for the fully 
amplifying (dot-dashed) and for lossy media (dotted curves) 
are similar as in Fig. [5] parameters for mixed-type (dashed 
curves) as in Fig. [l] The solid curve correspond to vacuum. 
The inset shows the corresponding a;-dispersion of e(iaj) = 
/i(ia;) for a mixed-type gain medium. 



In the main panel of Fig. [SJ showing the Casimir force 
as a function of distance, we see similar behavior: for 
low frequencies the curve is negative and the force is at- 
tractive but weaker than for free space, reminding of a 
purely lossy medium. For intermediate distances, the 
force is more attractive than for free space, as we have 
seen for fully amplifying media in Sees. IIVBI and IV Al 
Finally for large distances - but this is not clearly visible 
in the main graph - Fq again becomes weaker than i^™'^, 
again reminding of a lossy medium. 

These observations agree with the known fact that at 
long distances, the main contribution to the Casimir force 
comes from the low-frequency region, while the force 
at short distance depends on the high-frequency behav- 
ior [l^. This fact follows from the frequency- integral 
representation ([2^ of the force. From Fig. [3] for the fully 
amplifying media this relation was not evident, but in 
Fig. [5] it is: gain at a finite frequency interval around 
uJq may lead to Casimir forces that are more attractive 
than in free space for a finite interval of plate separa- 
tions, roughly around a separation d = Aq/ (27r), in terms 
of the resonance wavelength A§ = 2-kc/uJq. For closer or 
more distant separations, the lossy character of the gain- 
and-loss medium dominates the Casimir force, so that 
^vac < < 0. As for the fully amplifying media, we do 
not find Casimir repulsion for homogeneous media with 
both gain and loss. 



VI. CANONICAL QUANTIZATION OF 
ELECTROMAGNETIC FIELD IN AMPLIFYING 
MEDIUM 

In the previous sections we made a shortcut from the 
Lagrangian via the partition function and Green tensor 
to calculate the Casimir force for amplifying media. To 
find the Casimir force, we did not need to perform explic- 
itly a canonical quantization for amplifying media based 
on our Lagrangian. The reason to do this here is to make 
contact with other approaches and to provide an underly- 
ing canonical theory for the phenomenological quantum 
electrodynamics for amplifying media that was developed 
in recent years (33 - [38j |. 

The starting point is the Lagrangian ([1]), with the vec- 
tor potential A, and the continua of polarization oper- 
ators and Y„ as canonical fields with the following 
canonically conjugate fields 



eoE(x,i) 
n<,(x,i) 



5C 



(5A(x,t) 
5C 



eoA(x,t), 



(34a) 



(5X^(x,t) 

f{u, x) A(x, t) + sgn[£i(u;)]X^ (x, t), (34b) 

5C 
<5Y„(x,i) 

g(w,x)V X A(x,t) +sgn[^i(w)]Y^(x,i). 

(34c) 



So here we find — £oE as the canonical conjugate to 
the vector potential, as in Refs. [M-IH, l48l - [50| . Apart 
from the subtlety with the sign functions in Eqs. ([M)) 
that discriminate between the frequency intervals where 
there is gain and loss, the canonical quantization of the 
fields can proceed in a standard fashion by demanding 
equal-time commutation relations among the variables 
and their conjugates, 

[A,(x,0,-eo^j(x',0] =i^'5,j'5^(x-x'), (35a) 
[X^^,{-x,t),~'Q^, ^.j{-x! ,t)] = \hsg-D[£i{uj)]5ij 

X (5(w-cj')'5^(x-x'),(35b) 
[Y^^i{-K,t), -n^/j(x',t)] = \hsg-D.[m{uj)] 5ij 

X (5(a;-w')^^(x-x'),(35c) 

with all other equal-time commutators vanish. Using the 
Lagrangian ([1]) and the expression for the canonical con- 
jugate variables in p4p. we obtain the Hamiltonian den- 
sity 



n 



1 ^2/ N B2(x,t) 

-eoE2(x,t) ■ ^ ' ^ 



2^0 

) 

d^ Sgn[£i(L.)] {X2(x, t) 0.2X2 (x, 0} 

) 

Alo sgn[/.i {w)] {Yl (x, t) + uj^Yl (x, t) } . (36) 
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Maxwell's equations can now be obtained from the 
Heisenberg equations of motion for the vector potential 
and the transverse electric field and from the commuta- 
tion relation p5l) . 

(37a) 



A(x,i) = -E(x,t), 

• , , VxVxA(x,i) „ ^ v., 

eoE(x, t) = - V X M(x, t) - P(x, t). 



Mo 



(37b) 



Using the definitions D = EqE + P and H = B///o — M 
for the displacement field and the magnetic field strength, 
respectively, Eqs. (1571) result in D(x, t) = V x H(x, t) and 
B(x, t) — —V x E(x, t), as expected. In a similar fash- 
ion, the Heisenberg equation of motion for the dynamical 
variables X^^ and lead to 

X^(x,t) = -u?y.^{^,t)+^%^{ex{u:)\j{^,uj)E.{^,t),{Z%&) 
Y„(x,t) = -w2Y<,(x,i) -f sgn[^i(w)]5(x,w)B(x,i),(38b) 

with formal solution 

(x, t) = (x, 0) + X„ (x, 0) cos ujt 

+ /(x,c.)sgn[ei(^)] /'*dt'^^^^^^^^^^E(x,i'),(39) 
Jo ^ 

and likewise for Y^(x, t). To facilitate the calculations, 
let us introduce the following annihilation operators 



dj(x,a;,t) 
bj{x,u},t) 



1 



1 



[cjX^j(x,t)+iQ^,,(x,i)],(40a) 
[LjY^j(x,t)+in^j(x,t)],(40b) 



where j = 1,2,3 labels three orthogonal spatial direc- 
tions. Their commutation relations follow immediately 
from Eq. ([55]) . 



dj(x, w,t),d],(x',a;',t)l = sgn[£i(w)] 

X 5(cj- w')(5(x-x'),(41a) 

6j(x,w,t),6],(x',w',t) = sgn[^i(w)]%' 

X ,5(cj-w')'5(x-x').(41b) 

Now by inverting the relations (HUl) and substituting the 
result into Eqs. ([3]), the polarization and magnetization 
fields of the magnetodielectric medium can be written in 
terms of creation and annihilation operators as 

/>oo 

P(x, t)^eo dt' xc(x, t - t')E(x, t') + pN(x, i), 
Jo 

(42a) 

1 1"°° 

M(x, t) = — / di'xm(x, t - t')B(x, t') + MN(x, t), 
Mo Jo 

(42b) 



with susceptibilities Xm.c as defined in Eq. (|T5|) . The 
fields P^(x, t) and M^(x, t) are the electric and magnetic 



polarization noise densities associated with absorption 
and amplification. As in the phenomenological method, 
we can separate the noise operators into positive- and 
negative-frequency parts — P'^(+) -f pN(-) with 
pN(-) ^ [pN(+)]t ^j^jj analogously for M^, where 



/•OO 

P^(+)(x,0 = / dc^ 

"^O 



£oki(^^)| 



{d,(x,c^,0)e(ei(L^)) 
-f4(x,w,0)e(-£i(w))}e^''^*, (43a) 

M^(+)(x,t) ^ rd..W^^^M|6^(x,.;,0)e(^i(..)) 
Jo V '^Aio 



-h6t(x,c.,0)e(-w(c.))}£ 



(43b) 



In fact, the above equations are the starting point in 
Refs. (Sd- 38] to the phenomenological quantization of the 
electromagnetic field in amplifying magnetodielectric me- 
dia. 

If we now take the time derivative of Eq. (j37p and use 
Eq. (|42p . this yields the frequency-domain wave equation 
for the positive-frequency part of the vector potential 



V X [^"i(x,w)V X A(+)(x,a;)] - !^e(x, w)A(+'(x, 
-i^o^P^^+^(x) +MoV X MNW(x,tj). 
This equation can be solved as 

A(+)(x,t) = — / Alo I d3x'G(x,x',cj) • 
27r Jo J 

"-icjpN(+)(x',w)+V X MN(+)(x',w)" 



(44) 



(45) 



where the Green tensor G(x, x',aj) is the solution of 
Eq. (dH). 

The equations (|43|) and (|45|) and the commutation re- 
lations are the same as obtained from the phenomenolog- 
ical method (33 - [38l |. Therefore, with our Lagrangian ([T]) 
and the canonical quantization performed here, we for- 
mulated a microscopic basis for the phenomenological 
quantization of the electromagnetic field in amplifying 
magnetodielectric media. 



VII. CONCLUSIONS AND DISCUSSION 

The electromagnetic field in an amplifying magnetodi- 
electric medium was quantized with a path-integral tech- 
nique. We determined correlation functions of different 
fields and found electric and magnetic susceptibilities of 
the amplifying medium that are consistent with causality. 
We determined the Green functions in amplifying planar 
multilayer magnetodielectrics and used this to calculate 
the Casimir energy and force in such media. 

The calculations show that the form of the Casimir 
force as a functional of the dielectric functions 
£(x, w), /i(x, oj) does not change significantly as com- 
pared to passive media, but some more caution is needed, 
especially about the sign of wave-vectors components in 
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amplifying negative-index materials in the direction nor- 
mal to the planes. 

Here we studied systems that have both gain and a 
negative index, and to that end we assumed = 
in our numerical calculations, not because this would be 
simple to realize in experiments but rather to gain in- 
sight. The concept was introduced of fully amplifying 
media, i.e. amplifying at all frequencies. The advantage 
of our approach is that our optical functions e{uj),fi{uj) 
are defined at all frequencies and are causal. For two 
conductors separated by a homogeneous passive medium, 
we find that the Casimir force satisfies the bounds of 
Eq. (|77|) . so that it is always attractive and less so than 
if the medium were replaced by vacuum. On the other 
hand, for a homogeneous fully amplifying medium, we 
find the bounds (|26|l that the Casimir force is finite, and 
always more attractive than in vacuum. 

Both bounds are remarkable insofar that both the fully 
amplifying and the passive media with sicj) = may 
have negative refraction in a large frequency interval, 
whereas for some other planar geometries it was reported 
that negative refraction may lead to Casimir forces that 
may become repulsive [32[ , and also that gain would lead 
to repulsive Casimir forces [3l|, 111]. This is not a con- 
tradiction, however, but rather shows that negative re- 
fractive indices or gain do not automatically imply a sign 
change of the Casimir force, and that the sign of the force 
strongly depends on the geometry also for amplifying and 
negative- index materials. 

The bounds (|26l27p for homogeneous fully amplifying 
and for passive media have in common that the force is 
bound by the free-space Casimir force on the one side 
and by the free-space Casimir force divided by the static 
refractive index rigtatic on the other. It does not mat- 
ter whether for microwaves or optical frequencies the 
medium has a negative refractive index or not, as long as 
the static refractive index stays the same. In that sense, 
the Casimir force has little to do with optics, both for 
passive and for amplifying media. Ref. [25| also stresses 
the importance of the low-frequency behavior, and here 
we found an illustrative example also for amplifying sys- 
tems. 

Some observables have divergent values in models of 
linear amplification, especially for geometries with round- 
trip gain. Our simple geometry of perfect conductors sep- 
arated by a single amplifying medium will even exhibit 
round-trip gain for all frequencies for which there is am- 
plification. Nevertheless we find that the Casimir force 
on the conductors is finite, bounded by the inequalities 
of Eq. This could be explained by the fact that 

the Casimir force is a vacuum force, so that there are no 
photons present that are amplified indefinitely. 

Casimir forces in amplifying media are only beginning 
to be explored, and we only considered linear amplifi- 
cation. This is not always a realistic model, especially 
for geometries where the linear-amplification model pre- 
dicts round-trip gain. It is an open challenge to calculate 
Casimir forces in the presence of media with nonlinear 



amplification, for example with parametric amplification. 
A possible route could be to generalize the Lagrangian for 
passive nonlinear media of Ref. [67l | to gain media, and 
the result would be a generalization of our Lagrangian of 
Eq. ID) to nonlinearly amplifying media. 

Finally, we carried out a canonical quantization of the 
electromagnetic field in an arbitrary linear amplifying 
and/or passive medium, and showed that the resulting 
field operators satisfy the macroscopic Maxwell equa- 
tions for an arbitrary linearly amplifying and/or passive 
medium. The resulting theory is valid for all linear, inho- 
mogeneous, amplifying and/or passive magnetodielectric 
media with dielectric functions that satisfy the Kramers- 
Kronig relations. The postulates of the phenomenological 
theory [s^ - fssl ] that serve as its starting point were here 
derived by canonical quantization. 
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Appendix A: Generating function in presence of 
Dirichlet and Neumann boundary conditions 

Here we calculate the generating function for two per- 
fectly conducting plates surrounding an amplifying pla- 
nar multilayer system. We will be brief, as most of the 
calculation is identical to the case of lossy media, see 
Ref. ^ for example. 

We can consider TE and TM polarized waves sepa- 
rately. On the two plates (labeled by a) they satisfy 
Dirichlet or Neumann boundary conditions, respectively. 



fTM\so, =0, a = 1, 2 

Sn-^TElSo =0; a = 1,2 



(Ala) 
(Alb) 



on each surface Sa, where 9„ is the normal derivative 
of the surface Sa pointing into the space between the 
two plates. To obtain the partition function in 4D Eu- 
clidean space from the Lagrangian ([T|), we made a Wick 
rotation so that the signature of space-time changes from 
Minkowski to Euclidean. 

In the 4D Euclidean space, the plates are parameter- 
ized by A'i(x, zi,it) and A2(x, Z2,it), where x = {x,y). 
The Dirichlet or Neumann boundary conditions corre- 
sponding to the constraints (jAip can be imposed by in- 
serting delta functions which can be expressed in terms 
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of auxiliary fields ipaiXa) as follows |46l l68l| 



Using Eqs. (jA2|) . the partition function ([5]) in Euclidean 
space can be written as 

Ztm = Z^^ f V[^]l[S{^{X^))e-''^'™^'^\ (A3a) 

a=l 

where the Euclidean actions S'e,tm/te(v) S'I'g defined as 



+ /"dA'(5(A'-A'a)i/>„(a;)}, (A4a) 

'S'e.teM = j d'^x{C{ip{x)) 

+ (^(x) / dA'<5(A' - A'„)9„Va(a;)}. (A4b) 

a=l 

By comparing Eqs. (jA4p and (fTTj) . we can rewrite 
Eqs. dSS) as 

a = l 

^z(^J dX6{X - Xo,)ipa (A") j , (A5a) 

r ^ 

xZ dXS{X - A'„)a„?/'a('^)j ,(A5b) 

where the Z(. . .) in the two integrands are the generating 
functionals of interacting fields defined in Eq. ([TT]) with 
imaginary time. From Eqs. (jA5l) and (jl3l) the respective 
partition functions can be written as 



^TM/TE= / n^[V'o(^a)]e-''""'™^™''^°', (A6) 

where the effective actions SoS tm and SoS te are given 
by 

S'cff,TM(V'l,'02) = 

IY. 1 1 dX^dXp^^{X^)g{X^,Xfi)^p{Xp), (A7a) 



S'off,TE(l/'l7'02) 



\Y.jj dX^dXpM^o.)[dn^dn,g{X^,Xp)]^PpiXp). 

(A7b) 



Here the Green function of the fields after Wick rotation 
is denoted by the new font Q. The partition functions 
defined by (|A6|) are calculated straightforwardly, and the 
results are given in Eqs. ([^U)) and (PT|) in the main text. 



Appendix B: Green tensor for planar multilayer 
magnetodielectric media with gain 

For planar multilayer geometries as illustrated in 
Fig. m the electric permittivity and magnetic permeabil- 
ity vary only in the z direction, so we may introduce a 
transverse spatial Fourier transform as 

Gem(x - x', z, z', iuj) = J d^qe"''^''"'' ^GemIq, z, z' , iuj) 

(Bl) 

where q is a vector parallel to the conductor. Tomas 
uses this to arrive at the solution of Eq. in lossy 
dielectric multilayers [s^. The generalization to l ossy 
magnetodielectric media can be found in Refs. [sl, l69| . 
Here we briefly describe the results of a further nontrivial 
generalization, namely to magnetodielectrics with both 
loss and gain. In our notation we follow Ref. [56j . 

The Green tensor GEM(q7 z, z', iuj) assumes two differ- 
ent forms, depending on whether z and z' are located in 
the same layer or not. For z' in layer j it is given by 
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GEM(q, = 



GEM(q, z',iuj) 



jS{z — z')zz + 



TM 



-Qjdj 



=TE 3 



eo£j{\ijj)uj'^ 2Qj 

[£j^(q, iw; z)£j<(-q, iw; z')9(z - z') + £j<(q, iw; z)£j > (-q, iw; z')8(z' - z)], z in layer j 

(B2a) 

,. ^ TM ^o- „-(Qjdj + Q,d,) 



22/ 



£f >(q, ia;;z) „ 



£^<(-q,ic.;z')e(/-j) 



£f <(q,ia;;z) 



£^>(-q,ic.;z')e(j-0 



I/J 



, z in layer I ^ j 
(B2b) 



where ^te = —1, Ctm = 1, and Q{z) is the usual unit 
step function and 

£j>(q,iu;;z) - e+ (q)e-S.(--'^.) + rJ+e;^.(q)e^^-(^-'^\ 

(B3a) 
(B3b) 



r^_e+(q)e-e.-. 



£j<(q,ic.;z) = e-(q)ee.- 

Here a stands for TE or TM, and e^^^- = (q x z)j and 
^TMj ~ '^{^I'il^^ Qj^)j a-re the polarization vectors for 
TE and TM polarized waves propagating in the positive- 
/ negative-z direction, with qj = ^ w^Cj (ii^) / and 



for TE-polarized light, and for TM polarization 



„TM 



( ^J±l _ l2_\ 4- ( '^J±l I 



-2Qj±idj±i TM 



i±l± 



(B6b) 

For a finite number of layers there is only a finite number 
of relations to be solved, since for the leftmost and right- 

,1 1, ij J- 1 TE,TM n TE,TM n 

most layers one should take = 0, r„_|_ = 0, 

and di = dn = From the definition of the Fresnel 
coefficients introduced it follows that they satisfy 



,-2Qj±i<ij±i„TM 



Qj(q,iw) 



,(i^)M,(i^)/c2, (B4) 



which can be expressed in terms of the magnitude of the 
z-component Kj(q, cj) — ^u:'^ej(LS)^Xj(LS)j c? — cp- of the 
wave vector in layer j as Qj(q, iw) = — iKj(iaj). 

Here we arrive at a subtlety in the determination of 
the Green tensor for active multilayer media: the z- 
component of the wave vector, Kj (q, w), is not always well 
defined for real frequencies. The subtlety is that although 
the refractive index has no branch points in the upper 
half-plane, Kj(q, a;) may have branch points there 58]. 
If there are such branch points, then Kj(q, w) looses its 
usual physical interpretation. We will follow Refs. (58l[59| 
and only consider active media without branch points 



1 



'»/i/fc 



(B7) 



in the notation of Ref . [56| . 

We have hereby specified the rather complicated ex- 
pression for the Green tensor G of Eq. (jB2p . and we still 
need to relate it to the Green tensors G'^^'^^iq., z, z\ iuj) 
in terms of which the Casimir force is expressed in 
Eq. (j2ip . By using the ordinary coordinates according 
to the convention of Schwinger et al. Q and choosing 
(q, 2, q X z) — > (i, z, —y) the transverse electric and 
transverse magnetic Green functions satisfy 



d 



1 



d 



uj'^eiiijj, z) 



where is meaningful for real frequencies. In that 



dz n{\0J, z) dz 

case the signs of Ke[Kj{uj)] and Im[«;j(a;)] are identical to = /io5(z — z'), 

d 1 d 



those of Re[nj(a;)] and Im[nj(a;)], respectively, where rij 
is refractive index of j'-th layer (see Refs. [5^,|53|). Other 
quantities in Eqs. (IB2p that still need to be defined are 



1 



-2Q,di 



(B5a) 
(B5b) 



dz e(ia;, z) dz 
Hod{z - z'). 



fi{iuj,z) 



s{iuj, z) 



a;^/i(icLi, z) 



]g'^^{(i,z,z',iuj), 

(B8a) 

]c;™(q,z,z',ic.) 

(B8b) 



where rj_ and are the generalized coefficients for re- 
flection at the left / right boundary of layer j , which can be 
calculated with the aid of the recursive relations 63, £^ 



We checked with some lengthy but straightforward cal- 
culations that these Green functions ^'^^''^^(q, z, z', icj) 
with z' and z in layers j and I respectively, can be written 
very elegantly in terms of G(q, z, z', iw) as 



^TE 



( A'3±i _ tjS\ J_ ( W±i I h \ „-2Qj±i<ij±i TE 

VQj±i Qj ) ^ \ Qj±i ^ Qj J 

( ^j±i _L J_ f ^J±i _ tj.'] p-2Qj±i<ij±i„TE 

\Qj±i ^ Qj J ^ {Qj±i qJ^ ^J±i± 

(B6a) 



GTO(q,z,z',iw) = ^ (q,z,z',ia;) 

_ , , . , J. S(z — z') 
G^^(q, z,z ,10;) d,- 



e{iuj, z)£(ia;, z')a;^ 



£oe(iw, z)w2 
e™(q,z,z',ia.) 



(B9a) 



(B9b) 
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It is important to point out what has been achieved here: the equations (jB8l) , with boundary conditions that fohow 
by the identifications (IB9p we have found solutions for from the continuity of H^, Hy, and /xTJ^, and of E^, Ey, 
the scalar Green functions CJTE^tm g^j-g defined by and eE^- 
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